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Abstract 

The Bruggeman formalism provides an estimate e^J„ of the relative permittivity of a homogenized 
composite material (HCM), arising from two component materials with relative permittivities and eb- It 
can be inverted to provide an estimate of ea, from a knowledge of e£!„ and Numerical studies show that 
the inverse Bruggeman estimate ea can be physically implausible when (i) Re {e^Jto} /Rejcb} > and the 
degree of HCM dissipation is moderate or greater; or (ii) Reje^'j'^} / Re{eb} < regardless of the degree 
of HCM dissipation. Furthermore, even when the inverse Bruggeman estimate is not obviously implausible, 
huge discrepancies can exist between this estimate and the corresponding estimate provided by the inverse 
Maxwell Garnett formalism. 
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1 Introduction 

A composite material may be regarded as being effectively homogeneous, provided that wavelengths are 
much larger than the particle sizes of the component materials that make up the composite material. The 
constitutive parameters of such a homogenized composite material (HCM) can be estimated from a knowl- 
edge of the constitutive parameters of its component materials, along with a knowledge of the distributional 
statistics and shapes of its component particles [1, 2]. The Bruggeman homogenization formalism has been 
widely applied for this purpose for the past 70 years [1, 3]; and new areas of application for the Bruggeman 
formalism continue to emerge, for examples, in recent developments pertaining to complex HCMs [4, 5] and 
negatively-refracting metamaterials [2, 6]. However, a certain limitation of the Bruggeman homogenization 
formalism came to light in 2004 [7]. In the context of isotropic dielectric HCMs, derived from two compo- 
nent materials with relative permittivities ta and eh, it transpires that the Bruggeman estimate of the HCM 
relative permittivity may be physically implausible if ea/fb < in the case of nondissipative HCMs or if 
Re {ea} /Re {eb} < in the case of weakly dissipative HCMs. An example of such a problematic homogeniza- 
tion scenario — of considerable interest to the metamatcrial community — arises in the homogenization of 
silver particles with insulating particles at visible and near infrared wavelengths [8] . This limitation — which 
is relevant to active [9] as well as nondissipative and dissipative HCMs — also extends to the Maxwell Garnett 
homogenization formalism which shares a common provenance with the Bruggeman formalism [10], as well 
as the Hashin-Shtrikman, Wiener and Bergman-Milton bounds on the HCM's relative permittivity [7, 11]. 
Similar anomalous results, for HCMs arising from two component materials with Re{ea} /Re{e(,} < 0, have 
been described as 'electrostatic resonances' [12, 13, 14], but this term is avoided in [7, 8, 9, 11] (and herein) 
since the estimates of the HCM's relative permittivity described in [7, 8, 9, 11] are not physically plausible. 

Restricting our attention to the simplest possible case of an isotropic dielectric HCM arising from two 
isotropic dielectric component materials, in this communication we investigate the applicability of the Brugge- 
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man formalism to the inverse homogenization scenario wherein the relative permittivity of one of the compo- 
nent materials is estimated from a knowledge the relative permittivities of the other component material and 
the HCM. Formal expressions have been established for the inverse Bruggcman formalism (and the inverse 
Maxwell Garnett formalism) in the general setting of bianisotropic HCMs [15], but in certain cases these 
formal expressions may be ill-defined [16] and the ranges of applicability of these inverse formalisms have not 
been established. Our study is partly motivated by very recent implementations of the inverse Bruggeman 
formalism in estimating nanoscale constitutive and morphological parameters of certain sculptured thin films 
[17], which is a key step in modelling the electromagnetic response of infiltrated sculptured thin films [18, 19]. 



2 Analysis and numerical studies 

We consider the homogenization of two isotropic dielectric component materials with relative permittivities 
Ca and efc. The component materials a and h are assumed to be distributed randomly as spherical particles 
with volume fractions fa and fb = ^ — fa, respectively. The Bruggeman estimate of the relative permittivity 
of the corresponding HCM, namely e£r,„, is provided via [3] 

f - f^"^ fi. - f^-" 

f ^hcm I f ^hcm n /i ^ 

ta -r ^t/icm -r ^fc/jcm 

which is nonlinear in e^T^. A straightforward manipulation of (1) delivers the explicit formula 



{fa - 2/6) et + 2e£!;„ 

^hcm [■^) 



h {^b - + fa {e, + 2ef;„) 



for ea in terms of eb, e^J,„, fa and /;,. Since the component materials a and h are treated in an identical 
manner within the Bruggeman formalism, the corresponding formula for Cb has the same form as (2). Notice 
that as the inverse Bruggeman equation (2) does not involve a square root, there is no scope for Imjco} 
being nonzero if ef,, e^T^ 6 K- This contrasts with the forward Bruggeman formalism where a square root 
term enables Im {e£!m} to be nonzero even though e^, Cb G M. This physically-implausible scenario can arise 
when ea/eb < [7]. 

For comparison, we introduce the Maxwell Garnett estimate of the HCM relative permittivity [3] 



MG ^ , I 3/a£b (eg ~ £b) 

^hcm I ^ J. , 

Ea + 2eb - fa (ea ~ eb ) 



, 'JJa'^by-^a ^bl 

H -1^ n T y-^) 



and its corresponding inverse 



{2 + fa)e^,g~2fbeb 
{l + 2fa)eb-fbe^ig 



(4) 



The limiting behaviour of the inverse Bruggeman estimate (2) as compared with that of the inverse Maxwell 
Garnett estimate (4) is especially revealing. In the limit /a — > 1, both estimates yield the relative permittivity 
of the HCM, as they must.^ In the limit fa — >■ 0, the inverse Bruggeman formalism yields ea ~2efj^ 
whereas the inverse Maxwell Garnett formalism yields ea — 2eb. Therefore, the two inverse estimates differ 
markedly as fa approaches zero, provided that Cb and the relative permittivity of the HCM are sufficiently 
different. 

We now explore the inverse Bruggeman estimate (2), in comparison with the inverse Maxwell Garnett 
estimate (4), by means of some illustrative numerical examples. For nondissipative scenarios, the forward 
Bruggeman formalism runs into difficulties when ea/eb < but not when ea/eb > [7]. Accordingly, let us 



^Thc Maxwell Garnett estimate of the HCM relative permittivity is only strictly applicable in the dilute composite regime 
fa < 0.3. Accordingly, estimates of delivered by the inverse Maxwell Garnett formalism are strictly valid only for fa < 0.3. 
However, e^^J^ coincides with one of the Hashin-Shtrikman bounds on the HCM relative permittivity which applies at all values 

of fa [20], 
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begin by focussing on the regimes e^J,'„ /ef, < and e^J^ /ef, > 0. In Fig. 1, plots of Ca, as determined 
by the inverse Bruggeman formalism and the inverse Maxwell Garnett formalism, versus fa are provided 
for the cases where ei, = ±2 and — 3. When ^^cm^'^ /^b > the inverse Bruggeman and inverse 

Maxwell Garnett estimates arc in fairly close agreement. However, the values of yielded by the two 
inverse formalisms differ markedly when ^hcm^^ /^b < 0, except in the limit as fa approaches unity. Most 
notably, the inverse Bruggeman estimate becomes singular and undergoes a change in sign as the volume 
fraction increases through fa ~ 0.56, whereas the inverse Maxwell Garnett value remains finite and does not 
change sign. 

Next we turn to dissipative homogenization scenarios. In the case of the forward Bruggeman for- 
malism, problems arise when Re {cq} /Re {eb} < and the degree of dissipation is relatively small; if 
Re{ea} /Re {eh} < and the degree of dissipation is relatively large or if Rejca} /Re{e;,} > then the 
forward Bruggeman formalism was found to deliver physically plausible estimates of the HGM relative per- 
mittivity [7]. Accordingly, we consider the regimes where Re|e^J^^''| /Re{e6} > with the degree of 

dissipation in the HCM being relatively small, moderate and large. Graphs of the real and imaginary parts 
of Eq, as estimated by the inverse Bruggeman and inverse Maxwell Garnett formalisms, are plotted versus 
fa in Fig. 2 for the cases = 2 and e^^;^'^ = 3 + Si where S G {0.1, 1, 10}. When the de gree of HCM 
dissipation is relatively small {S = 0.1), the estimates of the real and imaginary parts of Ca provided by the 
inverse Bruggeman and inverse Maxwell Garnett formalisms agree fairly closely. When the degree of HCM 
dissipation is moderate {6 ~ 1), there is still fairly close agreement between the inverse Bruggeman and 
inverse Maxwell Garnett values of ea for most values of fa- Crucially, however, for fa < 0.05 the imaginary 
part of Ca estimated by the inverse Bruggeman formalism is negative-valued (unlike Im {ca} estimated by the 
inverse Maxwell Garnett formalism which is positive-valued). Here Im{ea} < is not a physically plausible 
outcome as it implies that the homogenization of an active material a and a nondissipative material b results 
in a dissipative HCM. For both the real and imaginary parts of the discrepancies between the values 
estimated by the two inverse formalisms become enormous when the degree of HCM dissipation is relatively 
large {S = 10). Furthermore, the inverse Bruggeman estimate is physically implausible for a much larger 
range of fa values; i.e., Im{ea} estimated by inverse Bruggeman formalism is negative-valued for fa < 0.3 
when 6 = 10. 

Lastly, we explore the Re |e^c'm^'^| /^^{^b} < regime. Plots of the real and imaginary values of 
in Fig. 3 correspond to the same parameter values as those used for Fig. 2 except that here et = —2. The 
estimates of the inverse Bruggeman formalism arc now physically implausible — due to Im {sa} < — for a 
wide range of fa values, regardless of whether the degree of HCM dissipation is relatively small, moderate or 
large. In contrast, the estimate of Im{ea} provided by the inverse Maxwell Garnett formalism is positive- 
valued for all scenarios considered. Additionally, the real parts of delivered by the two inverse formalisms 
differ enormously except when fa approaches unity, for all degrees of HCM dissipation considered. 

3 Closing remarks 

In the case of dissipative HCMs, the inverse Bruggeman estimates of can be physically implausible when: 

(i) Re {e£r,„} /Re {cb} > and the degree of HCM dissipation is moderate or greater; or 

(ii) Re {e£r„j } / Re {eb} < regardless of the degree of HCM dissipation. 

In the case of nondissipative HCMs. enormous discrepancies can exist between the estimates of ea provided 
by the inverse Bruggeman formalism and the inverse Maxwell Garnett formalism when ^hcm^'^ I ^b < 0. 
Therefore, the inverse Bruggeman formalism should be applied with great caution. Finally, we note that in 
the very recent implementations of the inverse Bruggeman formalism which motivated this study [17, 18, 
19], the relative permittivity parameters were positive-valued and the materials were nondissipative. The 
estimates yielded by the inverse Bruggeman formalism in these cases seem physically plausible, but the acid 
test can only be provided by suitable experimental measurements. 
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Figure 1: Plots of ta as determined by the inverse Bruggeman formalism (red, solid curves) and the inverse 
Maxwell Garnett formalism (blue, dashed curves) versus fa for = ±2 and ^^cm^^ — 3- Estimates of Ca 
delivered by the inverse Maxwell Garnett formalism are strictly valid only for fa ^ 0.3. 
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Figure 2: Plots of the real and imaginary parts of as determined by the inverse Bruggeman formalism 
(red, solid curves) and the inverse Maxwell Garnett formalism (blue, dashed curves) versus /a for = 2 and 
^hcm^^ = 3 + where 5 € {, 0.1, 1, 10}. Estimates of £a delivered by the inverse Maxwell Garnett formalism 
are strictly valid only for fa < 0.3. 
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